c 



Detecting the quantum zero-point motion of 
vortices in the cuprate superconductors 

Lorenz Bartosch Leon Balents*^, and Subir Sachdev'^ 

'^Department of Physics, Harvard University, Cambridge MA 02138 

^Institut fiir Theoretische Physik, Universitdt Frankfurt, 60054 Frankfurt, 

Germany 

Department of Physics, University of California, Santa Barbara, CA 93106-4030 



Abstract 

We explore the experimental implications of a recent theory of the quantum dy- 
namics of vortices in two-dimensional superfluids proximate to Mott insulators. The 
theory predicts modulations in the local density of states in the regions over which 
the vortices execute their quantum zero point motion. We use the spatial extent 
of such modulations in scanning tunnelling microscopy measurements (Hoffman et 
al, Science 295, 466 (2002)) on the vortex lattice of Bi2Sr2CaCu208+5 to estimate 
the inertial mass of a point vortex. We discuss other, more direct, experimental 
signatures of the vortex dynamics. 



1 Introduction 



It is now widely accepted that superconductivity in the cuprates is described, 
as in the standard Bardeen-Cooper-Schrieffer (BCS) theory, by the conden- 
sation of charge — 2e Cooper pairs of electrons. However, it has also been 
apparent that vortices in the superconducting state are not particularly well 
described by BCS theory. While elementary vortices do carry the BCS flux 
quantum of hc/2e, the local electronic density of states in the vortex core, 
as measured by scanning tunnelling microscopy (STM) experiments, has not 
been explained naturally in the BCS framework. Central to our considerations 
here are the remarkable STM measurements of Hoffman et al. [1] (see also 
Refs. [2,3]) who observed modulations in the local density of states (LDOS) 
with a period of approximately 4 lattice spacings in the vicinity of each vortex 
core of a vortex lattice in Bi2Sr2CaCu208+5. 

This paper shall present some of the physical implications of a recent theory 
of two-dimensional superfluids in the vicinity of a quantum phase transition 



Preprint submitted to Elsevier Science 



6 February 2008 



to a Mott insulator [4,5] (see also Ref. [6]). By 'Mott insulator' we mean 
here an incompressible state which is pinned to the underlying crystal lattice, 
with an energy gap to charged excitations. In the Mott insulator, the average 
number of electrons per unit cell of the crystal lattice, hmi, must be a rational 
number. If the Mott insulator is not 'fractionalized' and if umi is not an even 
integer, then the Mott insulator must also spontaneously break the space group 
symmetry of the crystal lattice so that the unit cell of the Mott insulator has 
an even integer number of electrons. There is evidence that the hole-doped 
cuprates are proximate to a Mott insulator with umi — 7/8 [7], and such an 
assumption will form the basis of our analysis of the STM experiments on 
Bi2Sr2CaCu208+5. The electron number density in the superfluid state, ns- 
need not equal hmi and will be assumed to take arbitrary real values, but not 
too far from Umi- 

A key ingredient in our analysis will be the result that the superfluid car- 
ries a subtle quantum order, which is distinct from Landau- Ginzburg order of 
a Cooper pair condensate. In two dimensions, vortices are point-like excita- 
tions, and are therefore bona fide quasiparticle excitations of the superfluid. 
The quantum order is reflected in the wavefunction needed to describe the mo- 
tion of the vortex quasiparticle. For umi not an even integer, the low energy 
vortices appear in multiple degenerate flavors, and the space group symmetry 
of the underlying lattice is realized in a projective unitary representation that 
acts on this flavor space. Whenever a vortex is pinned (either individually due 
to impurities, or collectively in a vortex lattice), the space group symmetry 
is locally broken, and hence the vortex necessarily chooses a preferred orien- 
tation in its flavor space. As shown in Ref. [4], this implies the presence of 
modulations in the LDOS in the spatial region over which the vortex executes 
its quantum zero point motion [8]. The short-distance structure and period of 
the modulations is determined by that of the Mott insulator at density umi, 
while its long-distance envelope is a measure of the amplitude of the vortex 
wavefunction (see Fig. 1). Consequently, the size of the region where the mod- 
ulations are present is determined by the inertial mass of the vortex. Here we 
will show how these ideas can be made quantitatively precise, and use current 
experiments to obtain an estimate of the vortex mass, rriy. There have been a 
number of theoretical discussions of using ECS theory [9,10,11,12,13,14,15], 
and they lead to the order of magnitude estimate ~ rrieikp^)'^, where rrie 
is the electron mass, kp is the Fermi wavevector, and ^ is the BCS coherence 
length. 



2 Vortex equations of motion 

We begin with a very simple, minimal model computation of the vortex dy- 
namics, in which retardation, dissipation, and inter-layer Coulomb interactions 
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Fig. 1. Schematic of the modulations in the LDOS of a vortex lattice. The short 

distance modulations in each vortex halo are determined by the orientation of the 
vortex in flavor space, as discussed in Ref. [4]. The envelope of these modulations 
is |^'(rj)p where ^ is the wavefunction of the vortices, and its characteristics are 
computed in the present paper. 

will be neglected. This serves the purpose of exposing the basic physics. The 
following section will present a much more complete derivation, in which these 
effects will be re-instated, and the connection to the field theory analysis of 
Refs. [4] will also be made explicit. 

Consider a system of point vortices moving in a plane at positions r^, where j 
is a label identifying the vortices. We do not exphcitly identify the orientation 
of each vortex in flavor space, because we are interested here only in the long- 
distance envelope of the LDOS modulations; the flavor orientation does not 
affect the interactions between well-separated vortices, and so plays no role 
in determining the wavefunction of the vortex lattice. In a Galilean-invariant 
superfluid, the vortices move under the influence of the Magnus force 

where t is time, h — 27ih is Planck's constant, Vs(r) is the superfluid veloc- 
ity at the position r, and ns/o^ is the electron number density per unit area 
(a^ is the area of a unit cell of the underlying lattice). One point of view is 
that the force in Eq. (1) is that obtained from classical fluid mechanics after 
imposing the quantization of circulation of a vortex. However, Refs. [16] em- 
phasized the robust topological nature of the Magnus force and its connection 
to Berry phases, and noted that it applied not only to superfluids of bosons, 
but quite generally to superconductors of paired electrons. Here, we need the 
modiflcation of Eq. (1) by the periodic crystal potential and the proximate 
Mott insulator. This was imphcit in the results of Ref. [4], and we present it 



in more physical terms. It is useful to first rewrite Eq. (1) as 



m„^ = F£(j) + FB(j), (2) 

where is the first term proportional to and F^ is the second term. Our 
notation here is suggestive of a dual formulation of the theory in which the 
vortices appear as 'charges', and these forces arc identified as the 'electrical' 
and 'magnetic' components. In the Galilean invariant superfiuid, the values 
of F^' and F^ are tied to each other by a Galilean transformation. However, 
with a periodic crystal potential, this constraint no longer applies, and their 
values renormalize differently as we now discuss. 

The influence of the crystal potential on F^; is simple, and replaces the number 
density of electrons, ns, by the superfiuid density. Determining Vs(rj) as a sum 
of contributions from the other vortices, we obtain [17] 

FM - ^^Ps E , (3) 

where Ps is the superfiuid stiffness (in units of energy). It is related to the 
London penetration depth. A, by 

^^~167re2A2' 
where d is the interlaycr spacing. 

The modification of F^ is more subtle. This term states that the vortices are 
'charges' moving in a 'magnetic' field with ns/2 'fiux' quanta per unit cell 
of the periodic crystal potential. In other words, the vortex wavefunction is 
obtained by diagonalizing the Hofstadter Hamiltonian which describes motion 
of a charged particle in the presence of a magnetic field and a periodic po- 
tential. As argued in Ref. [4], it is useful to examine this motion in terms of 
the deviation from the rational 'fiux' Hmi/'^ = p/q (Pi 9 are relatively prime 
integers) associated with the proximate Mott insulator. The low energy states 
of the rational flux Hofstadter Hamiltonian have a g-fold degeneracy, and this 
constitutes the vortex flavor space noted earlier [18]. However, these vortex 
states describe particle motion in zero 'magnetic' field, and only the deficit 
[ns — Umi)/'^ acts as a 'magnetic fiux'. This result is contained in the action 
in Eq. (2.46) of Ref. [4] (see also Eq. (17) below), which shows that the dual 
gauge flux fluctuates about an average flux determined by {ng — umi)- The 
action in Ref. [4] has a 'relativistic' form appropriate to a system with equal 
numbers of vortices and anti-vortices. Here, we are interested in a system of 
vortices induced by an applied magnetic field, and can neglect anti- vortices; so 
we should work with the corresponding 'non-relativistic' version of Eq. (2.46) 
of Ref. [4]. In its first-quantized version, this 'non-relativistic' action for the 
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vortices leads to the 'Lorentz' force in Eq. (2) given by 



I ■^ K'^s - nui) drj , 
2^^ 

If the density of the superfluid equals the commensurate density of the Mott 
insulator, then = 0; however, remains non-zero because we can still 
have Ps 7^ in the superfluid. These distinct behaviors of Fe,b constitute a 
key difference from Galilean-invariant superfluids. In experimental studies of 
vortex motion in superconductors [19], a force of the form of Eq. (5) is usually 
quoted in terms of a 'Hall drag' co-efficient per unit length of the vortex line, 
a; Eq. (5) implies 

h{ns - umi) 

" = 2^^- 

Thus the periodic potential has significantly reduced the magnitude of a from 
the value nominally expected [20] by subtracting out the density of the Mott 
insulator. A smaller than expected |q;| is indeed observed in the cuprates [19]. 
It is worth emphasizing that F^ (but not Fg) is an intrinsic property of a 
single vortex. Moreover, we expect that, taken together, the relation Eq. (6) 
and the fiavor degeneracy q are robust "universal" measures of the quantum 
order of a clean superconductor, independent of details of the band structure, 
etc. 



3 Derivation from field theory 



Wc will now rederive the results of the previous section from a more sophisti- 
cated perspective. We will use a field theoretic approach to derive an effective 
action for the vortices, a limiting case of which will be equivalent to the equa- 
tions of motion already presented. The effective action will include retardation 
effects, and can be easily extended to include inter-layer interactions and dis- 
sipation. 

Our starting point is a model of ordinary bosons on the square lattice interact- 
ing via the long-range Coulomb interaction. Following Ref. [4] we will briefly 
review a duality mapping of this model into a field theory for vortices in a 
superfluid of bosons which is in the vicinity of a transition to a Mott insulator. 
The density of bosons per unit cell of the underlying lattice is ps = riB/a^, 
while the density of bosons in the Mott insulator is pmi — ''t-mI/o^ = {p/(l)/ct'^', 
here is the unit cell area of the underlying lattice. Closely related held the- 
ories apply to models of electrons on the square lattice appropriate to the 
cuprate superconductors [5,21], with the boson density replaced by the corre- 
sponding density of Cooper pairs; the needed extensions do not modify any of 
the results presented below. It should be noted that since two electrons pair to 
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form one Cooper pair the average number of electrons in the Mott insulating 
state is umi = '^^Umi (and the average number of electrons in the superfluid 
phase is Ug = 2ns). 

In zero applied magnetic field, the Hamiltonian of our system is given by 

u \^ fA 1\ < e*^ ^{ni-nB){nj -ub) 

W = -p.2^cos(Ac,0ij + — }^ 1 . _ ^ , (7) 

ia i^j I * 3\ 

where ps is the superfluid stiffness and — e* (= — 2e) is the charge of a boson 
(Cooper pair). The bosons are represented by conjugate rotor and number 
operators (pi and hi which live on the sites i of the square lattice (with position 
vector Ti) and satisfy the commutation relations 

[(i)i,hj\^i5ij . (8) 

We subtract the average boson density ub from the number operators fii to 
account for global charge neutrahty of the system. Finally, we have introduced 
the discrete lattice derivative A^^j = <t>i+a ~ (t>i along one of the two spatial 
directions a — xory. 



3.1 Dual lattice representation 



Let us now briefly review the duality analysis of the above model with spe- 
cial emphasis on the long-range Coulomb interaction. Following Ref. [4] we 

represent the partition function of ?i as a Feynman path integral by inserting 
complete sets of eigenstatcs to the number operators hi at times separated by 
the imaginary time slice Ar. While the Coulomb interaction term is diagonal 
in this basis, the hopping term in Ti. can be easily evaluated by making use of 
the Villain representation 

exp (p^Arcos (Aa0i)) ^ cxp [ - "^'^ iJi^l^^cpA . (9) 

Here, we have set = 1 and have dropped an unimportant normalization 
constant which we will also do in the following. The JjQ arc integer variables 
residing on the links of the direct lattice, representing the current of the bosons. 

Extending the lattice index i to spacetime and introducing the integer- valued 
boson current in spacetime, Jj^ = (rij, Jj-^, Jj^), the partition function can be 
written as 
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where the prime on the sTim over the Jj^ restricts this sum to configurations 
satisfying the continuity equation 



A^J,^ = . (11) 
This constraint can exphcitly be solved by writing 

Jifi = ^fiuX^u^ix , (12) 

where Aj^j^ is an integer- valued gauge field on the links of the dual lattice with 
lattice sites i. We can now promote Aj^ from an integer- valued field to a real 
field by the Poisson summation method. We then soften the integer constraint 
with a vortex fugacity yy and make the gauge invariance of the dual theory 
explicit by by replacing I'kAj^ by 2TiAi^ — A^t?^. The operator e^'^* is then the 
creation operator for a vortex in the boson phase variable 0^. We now arrive 
at the dual partition function 



^'i = ]JiJ dAi^ j d-di exp ^ cos (2?!^^^ - A^^?^) 



2Arp, 

Are*^ ^ {eo^xKA-jx - nB){eQ^'\'K'A-jx' - ub) i ^^^^ 



As a last step we can replace the hard-core vortex field e*''' by the "soft-spin" 
vortex field t/^, resulting in 



= n / dAi, j d^k exp ( I E [^U,e'^'^'^^k + cc] - ^ 



+ 



u 



2Atp, 



Are*^ ^ {eoi^x^vAix - nB){eoi,'x'^v'A-jX' - ^b) 



E 



r,: - r,- 



(14) 



The first two terms in the exponent describe the action of the vortex fields ■^j 
which are minimally coupled to the gauge field A^^. While the system is in a 
superfiuid phase for s it is in a Mott insulating phase for s <C 0. 

At boson filling n\^j = p/q the gauge field Aj^ in the action in Eq. (14) fiuc- 
tuates around the saddle point ^4^^ with Cnux^uAjx — '>^Mi^iJ.,r- It is therefore 
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customary to substitute the gauge field by 

Ar 

■^ir ~^ ^It + "2^^^"^ ' ^^^^ 
Ma ~^ Ma + 2^Ma ■ (16) 

Here we have aheady rescaled the deviations of the gauge field from the sad- 
dle point such that later on we can easily take the continuum limit. A care- 
ful analysis of the symmetry properties of the above dual vortex theory (see 
Ref. [4]) shows that the vortex fields transform under a projective symmetry 
group whose representation is at least g-fold degenerate. It was also argued in 
Ref. [4] that while q cannot be chosen too large the boson density in the su- 
perfluid phase ub can take any value not too far away from the boson density 
in the Mott insulating phase, Umi- 

In zero applied magnetic field, and at a generic boson density pb, the field 
theory for such a superfiuid is then given by 



=0 



g*2 



X 



(z • (V X A(r, t)) - 27r(pB - pui)) (z • (V x A(r^ r)) - 27r(pg - pmi)) 

|r — r'l 

+ ... (17) 

This equation is a modified version of Eq. (2.46) in Ref. [4] with the short-range 
interaction between bosons replaced by the long-range Coulomb interaction. 
ip£ is a vortex field operator which is the sum of a vortex annihilation and an 
anti-vortex creation operator, and i is the vortex flavor index. As discussed 
above, as long as the vortices are well separated, the flavor index i plays no role 
in determining the zero-point motion of the vortices, and hence the envelope 
of the modulations illustrated in Fig 1; we will therefore drop the flavor index 
in the subsequent discussion. Recall that the index p runs over the spacetime 
co-ordinates r, x, y (while the index a runs only over the spatial co-ordinates 
X, y). We have rescaled the r co-ordinate so that the 'relativistic velocity' 
appearing in the flrst term is unity. 

The vortices in are coupled to a non-compact U(l) gauge field ~ 
(A,-, A). The central property of boson- vortex duality is that the 'magnetic' 
flux in this gauge field, z • ( V x A)/ (27r) is a measure of the boson density. 
However, notice from the last term in S^p with co-efficient e*^ that the action 
is minimized by an average gauge flux (or boson density) of (p^ — Pmi), the 
deviation in the density from that of the Mott insulator, and not at the to- 
tal boson density pb, as one would expect from usual considerations of the 
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Magnus force on continuum superfluids. The origin of this shift in the average 
flux is explained in detail in Ref. [4]; briefly stated, the combination of the 
periodic potential and 'magnetic' flux acting on the vortices has the effect of 
transmuting the flux associated with the density of the Mott insulator into 
the q vortex flavors. Only the deflcit ps — Pmi then acts as a 'magnetic' held 
on the vortices. 



The vortices also experience an 'electric' fleld, whose fluctuations are con- 
trolled in the action by the boson superfluid density Ps- 



The vortex component of the action has a 'relativistic' form and so de- 
scribes both vortices and anti-vortices with vanishing net mean vorticity. We 
are interested here in the case of a vortex lattice induced by an applied (real, 
not dual) magnetic fleld. In the dual language, this magnetic fleld appears as a 
static background 'charge' density which interacts via the 'electric' force with 
the 'charged' vortices and ant i- vortices. Finiteness of energy requires that this 
background charge density induces a neutralizing density of vortex 'charges', 
which, in the classic Abrikosov theory, form a vortex lattice (in the dual lan- 
guage this lattice is a Wigner crystal of charges) . We will neglect anti- vortices 
from now on, and focus only on the dynamics of these vortices induced by the 
applied fleld. For the action S^p this restriction means that we should work 
with the 'non-relativistic' limit. The formal procedure for taking this limit 
was discussed in Section IV. B of Ref. [4], and leads to an action for a non- 
relativistic fleld which is a vortex annihilation operator (anti- vortices have 
been eliminated from the spectrum). As shown earlier, the action for ^ takes 
the form 



2m„ 



Stt^Ps 



+ 



,*2 



87r2 



I d\ I dh' j dr 



X 



(z • (V X A(r, r)) - 2t:{pb - Pmi)) (z • (V x A(r', r)) - 2t:{pb - Pmi)) 



r — r 



+ ... 



(18) 



We now transform from this second quantized form of the vortex action to a 
flrst quantized form with vortices at spatial positions rj(r) where, as before. 
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j is a vortex label. In this form the action is 



2 



+ 



g*2 



I (fr j £r' I dr 



X 



87r2 

(z • (V X A(r, t)) - 27r(pB - pmi)) (z • (V x A(r^ r)) - 27r(pg - pmj)) 

|r — r'l 

+ ... (19) 
where p and J are vortex density and currents respectively: 

= X!^(r- rj(r)) , 
J(r,T) = E^'^(r-r,(r)). (20) 

Now it is useful to shift the vector potential A^, to absorb the mean background 
flux 

A^^B + A, (21) 

where 

B = 2n{pB - pmi) ■ (22) 
The fluctuations of the flux about this average value are controlled by the 
long-range Coulomb interactions. We assume that the vortices are located 
near the positions of a regular vortex lattice with equilibrium positions Rj, 
and make displacements Uj from these positions such that = Rj + u^. 
Adopting the Coulomb gauge, V • A = 0, the resulting action for the vortices 
is Su = Si + S2 + S3 where 



duj \ ^ B ^ ( duj 

1^ +^2^-r^-"^ 



*2 

2 



S3 = J d^rdr {iA^p + iJ ■ A) . (23) 

It is interesting to note that all couplings in this action are known from ex- 
periments, apart from the vortex mass m^. 

Now we integrate out the Aj- and A, and expand the resulting action carefully 
to second order in the u. (We also use the component notation m^, where the 
index a extends over the x and y components.) This directly yields the result 



X p duo f d/'^ 
^-2^1 ^ J,^^^2M-CL,-u^)D^pMu,{<i,u;) , (24) 
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where the momentum integral is over the first Brillouin zone of the vortex 
lattice, 

u(q,u;) = J drY.^je''"^'^'^"'"^ , (25) 



and the dynamical matrix is 



Docp (q, ^^) = AomyUj'^docp + AoUjBeap , 

G^O 1^ 

^ |q + G|2 + cu2|q + G|/(27rp3e*2) 



G 

where Aq is the area of a unit cell of the vortex lattice, and G extends over 
all the reciprocal lattice vectors of the vortex lattice of points Rj. 

It is now not difficult to show (see Appendix A) that, after dropping retarda- 
tion effects, the action in Eqs. (24, 26) is equivalent to the harmonic equations 
of motion that would be obtained for the vortex lattice from Eqs. (2-5). In- 
stantaneous interactions are obtained by taking the e* — oo limit of Eq. (26). 
Clearly, the present formalism allows us to include these without much addi- 
tional effort. 

So far, the action is free from dissipation effects associated with the Bardeen- 
Stcphcn viscous drag. We will consider these in Section 5.1 below. For now we 
note that these can be included in the above action simply by the transfor- 
mation 

As we will see in Eq. (36) , rj is the viscous drag co-efficient, and d is the spacing 
between the layers. 

The present formalism also allows us to consider the coupling between different 
two-dimensional layers in the cuprate system, and this will be examined in the 
following subsection. 



3.2 Interlay er Coulomb interactions 



Even in the absence of any Josephson or magnetic couplings between the 
layers, it is clear that we at least have to account for the interlayer Coulomb 
interactions because the vortex spacing is much larger than the layer spacing 
d. 
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We create a copy of all fields in all layers, labelled by the layer index n. In 
particular, we now have gauge fields A^^\ The Coulomb couplings between 
the layers modify ^2 to 



cPq / 1 



27r 7iBZ 47r2 \ Src^p 



s n 



,*2 



+ ^Y1 e-l"-"'l«'^(q X A(")(-q, -co)) ■ (q x A("')(q, co)) ] . (28) 



The intcrlayer interaction comes from the Fourier transform of 1/y + (n — n')'^d?. 
Now we perform a Fourier transform of the layer index, into a momentum per- 
pendicular to the layer, p^, leading to the field >l^(q, a;,px). In terms of this 
field 



So 



^ dp± f duj f d'^q ( 1 



TT 27: J 27r iiBZ 47r^ I Sn'^ps „ 



^ [g (q, a;, px) I ' + c^' I A (q, a; , px) P 



,*2 



+ 



1 - e-2«<i 



Anq 1 + e-29d - 2e-«<^cosp_L 



(q X A(-q, -u, -p±)) ■ (q x A(q,w,px)) | . 

(29)^ 



Because we always have go? <S 1, we can simplify this to 



dp± f du f d^q I 1 



+ 



TT 2tt J 2n JiBZ Att^ \ Stt^p^ 



J2 [q^\A^iq,u;,p±)\'' + u;^\A{q,uj,pj_)f 



47r(l — cosp_l) 



q X A(-q, -a;, -p±)) ■ (q X A(q,a;,px)) 



(30) 



Now, as before, we integrate out the and obtain the effective action for the 
u in a single layer. This has the form as in Eq. (24), but with the dynamical 
matrix in Eq. (26) replaced by 



Da/siq, uj) = AomyUj'^Saf} + AouBea/s - X] i^^'^'^ 

dp± A'K^Ps{qa + Ga){qf3 + Gf3) 



+E 

G 



Q .-TT 27r |q+G|2 + u;2(l -cosp±)/(27rp,e*2d) 
dp± An'^psou^il — cosp±) 



G 



-TT 27r a;2(l -cospx) +27rpse*2|q + G|2d ' 



(31) 
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We begin by presenting the numerical solution of the minimal model presented 
in Section 2. This is equivalent to solving the dynamical matrix in Eq. (26) 
in the limit e* — oo. The influence of all the additional effects considered in 
Section 3 will be described in the next section. 



We evaluated the dynamical matrix for a perfect triangular lattice of vor- 
tices at positions using the Ewald summation technique (see Appendix 
A). This leads to the vortex 'phonon' modes shown in Fig. 2 and vortex 'mag- 
netophonon' modes shown in Fig. 3. The computation of these modes is a 
generalization of other vortex oscillation modes discussed previously in super- 
conductors [22,23], rotating superfluids [24,25,26], and, in a dual picture of 
'charges', also to oscillations of electronic Wigner crystals in a magnetic field 
[27,28]. 

Quantizing these modes, we determine the mean square displacement of each 
vortex due to the quantum zero point motion of the vortex lattice, which we 
denote u^^^ = (jr^ — Rjp)/2. In terms of the normal modes we find 

, h 

7/^ — 

2m„a;, 



+ 



ujrjAr, r d q 

X ^ ' 



2 Abz (27r)^ 



,<^i(q) <^2(q) a;i(q)a;2(q)[a;i(q) -Fa;2(q)] 



(32) 
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Fig. 3. Dispersion of the 'magnetophonon' modes of the vortex lattice (here we take 



0.5a;„ 



Here, the momentum integral is over the first Brillouin zone. The prefactor 
h/2myUip should be identified as the mean square deviation of the position 
of a simple one-dimensional oscillator of mass and frequency ujp from its 
equilibrium position. 



We found an excellent fit (see Fig. 4) of our numerical data to the interpolation 
formula 

0i)3627a|^ ^ 

m„ = 1 F{x), (33) 

where is the separation between nearest neighbor vortices, x = \a\du. 
and 

Fix) ^ 1 - 0.4099a;^ . 



2 
rms 



/h, 



(34) 

The above fit is motivated by a simple analytic calculation in which the two 
modes a;i,2(q) are replaced by their long- wavelength approximations (see Ap- 
pendix B for details). Eq. (34) holds only as long as the r.h.s. is positive, while 
F{x) — for larger x (we will see below and in Fig. 5 that this apparent upper 
bound on x is relaxed once we allow for viscous damping). Similarly, for the 
frequency cUmp in Fig- 3, we obtain cu^p = + uj^ with 



Up = 



2 

rms 



35.45psW 



21.^1psu'i^^xF{x) 



(35) 



For the experiments of Ref. [1] wc estimate Ps = 12 meV [29], Mrms = 20 A 
[30], a = 3.83 A, and = 240 A. The overall scale for rriy is determined by 
setting ns — umi so that x — and F{x) = 1. This yields rriy ~ Snie and 
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Fig. 4. Plot of our numerical data (open circles) and fit to our interpolation formula 
F{x) (straight line) as given by Eq. (34). 

Tujjp K, 3 meV (or z/p pa 0.7 THz). For a more accurate determination, we need 
nsi for which there is considerable uncertainty e.g. for \ns — umi] — 0.015, we 
find X — 1.29, m„ ~ 3me and liujp ~ 5 meV. 



5 Limitations 

We now consider the influence of a variety of effects which have been neglected 
in the computation of Section 4. The extensions were already discussed in 
Section 3, and here we will make quantitative estimates. 

5.1 Viscous drag 

It is conventional in models of vortex dynamics at low frequencies [31] to 
include a dissipative viscous drag term in the equations of motion, contributing 
an additional force 



to the r.h.s. of Eq. (2). This leads to the transformation Eq. (27) in the dy- 
namical matrix. There are no reliable theoretical estimates for the viscous 
drag co-efficient, ?], for the cuprates. However, we can obtain estimates of 
its value from measurements of the Hall angle. Oh-, which is given by [31,32] 



(36) 
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Fig. 5. Plot of the function F{x,y) which replaces F{x) in Eqs. (33,34,35) upon 
including viscous drag, rj (y = rjdu^^^^/fi). The argument x measures the Hall drag 
a {x= \a\du'^^^/h) and the Hall angle is determined by | tan6^| = x/y. 



tanOn — ck/tj. Harris et al. [32] observed a dramatic increase in the value 
|tan^/f|, to the value 0.85, at low T in "60 K" YBa2Cu306+j/ crystals, sug- 
gesting a small t], and weak dissipation in vortex motion. For our purposes, we 
need the value of r] for frequencies of order 0;^, and not just in the d.c. limit. 
The very naive expectation that 77(0;) behaves like the quasiparticle microwave 
conductivity would suggest it decreases rapidly beyond a few tens of GHz, well 
below Up [33] . Lacking solid information, we will be satisfied with an estimate 
of the influence of viscous drag obtained by neglecting the frequency depen- 
dence of rj (a probable overestimate of its influence). The resulting corrections 
to Eqs. (33-35) are easily obtained (as in Ref. [11]), and can be represented 
by the replacement 

F{x) F{x, y), where y = V^l^ (37) 

and 

F{x,Q)^F{x). (38) 

The sketch of the function F{x,y) is in Fig. 5; as long as y > 0, we have 
F{x,y) > 0. As expected, the viscous damping decreases the estimate of the 
mass, and this decrease is exponential for large y, e.^. at x = we have the 
interpolation formula 

F(0, y)^(l + 0.41y + 2.69?/^) e^^^'^^^. (39) 
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5.2 Meissner screening 



The interaction in Eq. (3) is screened at long distances by the supercurrents, 
and the intervortex couphng becomes exponentially small. This does have an 

important influence at small momenta in that the shear mode of the vortex 
lattice disperses as [22,23] ~ q^. However, as long as <C A, there will not be 
a significant infiuence on Wrms oi' ^p- 

5.3 Retardation 

The interaction Eq. (3) is assumed to be instantaneous; in reality it is retarded 
by the propagation of the charged plasmon mode of the superfluid, and these 
effects were included in Eqs. (26) and (31). We can estimate the corrections due 
to this mode in a model of superfluid layers coupled by the long-range Coulomb 
interaction. In physical terms, we compare the energy per unit area of a 'phase 
fluctuation' at the wavevector of the vortex lattice Brillouin zone boundary 
(~ Ps/al) with its electrostatic energy (~ h'^ujp/{e'^d)); this shows that such 
corrections are of relative order ~ {fi^ / m,^) / {e d) . Alternatively this ratio can 
be viewed as the order of magnitude of the two terms in the denominator of 
the fourth term in Eq. (31). For the parameters above and d = 7.5 A this ratio 
is ~ 0.009, and hence quite small. 

5.4 Nodal quasiparticles 

We expect that nodal quasiparticles contribute to the viscous drag, and so their 
contribution was already included in the experimentally determined estimate 
of r/ in (i). The nodal quasiparticle contribution to and f] has recently 
been discussed at some length in Ref. [34]. This analysis flnds an infrared 
flnite correction to and a contribution to rj which vanishes as T — > 0. The 
latter observations are consistent with the observations of Harris et al. [32] . 

5. 5 Disorder 

We have assumed here a triangular lattice of vortices. In reality, STM experi- 
ments show signiflcant deviations from such a structure, presumably because 
of an appreciable random pinning potential. This pinning potential will also 
modify the vortex oscillation frequencies and its mean square displacement. 
Both pinning and damping rj tend to reduce vortex motion. For this reason, the 
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estimates of m„ above in which these effects are neglected must be regarded 
as upper bounds. 

The above considerations make it clear that new experiments on cleaner un- 
derdoped samples, along with a determination of the spatial dependence of 
the hole density (to specify a), are necessary to obtain a more precise value 
for m„; determining the H dependence of rriy will enable confrontation with 
theory. 



6 Implications 

An important consequence of our theory is the emergence of cUp as a character- 
istic frequency of the vortex dynamics. It would therefore be valuable to have 
an inelastic scattering probe which can explore energy transfer on the scale of 
hijp, and with momentum transfer on the scale of h/uy, possibly by neutron 
[35] or X-ray scattering; observation of a resonance at such wavevectors and 
frequencies, along with its magnetic field dependence, could provide a direct 
signal of the quantum zero-point motion of the vortices. A direct theoretical 
consideration of magnetoconductivity in our picture would have implications 
for far-infrared or THz spectroscopy, allowing comparison to existing experi- 
ments [36] ; further such experiments on more underdoped samples would also 
be of interest. 

Another possibility is that the zero point motion of the vortices emerges in 
the spectrum of the LDOS measured by STM at an energy of order hujp. We 
speculate that understanding the 'vortex core states' observed in STM studies 
[2,3] will require accounting for the quantum zero point motion of the vortices; 
it is intriguing that the measured energy of these states is quite close to our 
estimates of hujp. 
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A Evaluation of 'magnetophonon' modes 



In this appendix we show how the 'phonon' and 'magnetophonon' modes as 
depicted in Figs. 2 and 3 can be calculated using the well-known Ewald sum- 
mation technique (see e.g. Ref. [37]). Our calculation is similar to existing 
calculations and we will make contact with earlier work by considering a gen- 
eralized potential energy of the form 

with p > 0. Here, is the position of the i'th 'charge' (which could also be a 
real charge) and denoting Rj as is its equilibrium position we write = Rj-t-Uj. 
To minimize the total potential energy the Rj form a triangular Bravais lat- 
tice. Forp = 1, [/ reduces to the potential energy of a two-dimensional Wigner 
crystal with charges interacting via the three-dimensional Coulomb interac- 
tion. We are particularly interested in the case p ^ where the interaction 
between 'charges' becomes logarithmic. This case applies to our vortex lat- 
tice: Taking the gradient of U with respect to rj, we obtain the 'electric' force 
F£;(j) = — Vr^C/ on the j'th vortex from 



after setting p — 0. Identifying g'^ — 27rps, this equation clearly reduces to 
Eq. (3). 

By considering arbitrary p we will now generalize a calculation by Bonsall and 
Maradudin [28]. First, we expand U in the displacements from the equilibrium 
positions (with a — x,y labeling the two cartesian coordinates) and only 
keep terms up to second order, 

U = Uo + ^Yl UiaUjf3 ■ (A.3) 

To determine the normal modes we essentially have to diagonalize the matrix 

Fourier transforming ^fa-j/s to momentum space gives us a block-diagonal ma- 
trix fi^^(q) where each block is a 2 x 2 matrix. The action for the vortices is 
then given by Eq. (24) with the dynamical matrix 

Da/3{q, uj) = Ao'iJiyOu'^Sap + AoBuea/3 + AorriyQlJci) . (A.5) 
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Due to the long-range interaction all matrix elements fi^^(q) are slowly con- 
verging sums which we evaluate using the Ewald summation technique [37]. 

Let us write r2^^(q) as 



fill) 



(A.6) 



where the matrix elements 5*0/3 (q) are defined as 



E 



tq-R, 



x=0 



(A.7) 



We can now use the integral representation 



1 



pyP 2r(l+p/2)7o 



(A.8) 



(with arbitrary e > 0) and divide the integral on the r.h.s. in one part with 
< t < 1 and one part with t > 1. Setting y = |x — Rj| we see that the sum 
over the Rj in the integral from 1 to infinity converges rapidly. The usual trick 
is to use Ewald's generalized theta function transformation 



- iq- Rj — I X — Rj p et 



TT 



Aoet 



E 

G 



,i(G+q)-x-|G+q|2/4et 



(A.9) 



and transform the integrand of the integral from to 1 to a sum over the 
reciprocal lattice G such that this sum also converges rapidly. We then obtain 



5a/3(q) 



21(1 +P/2) 



+ 



Lr 

4e5„ a IT 



2+p Ao€ 



E(?« + G^){qp + Gf3)ip_p/2{\q + G|74e) 



G 



(A.IO) 



where 



dtfe 



V —zt 



(A.11) 



20 



is a Misra function. The matrix Jla^(q) is now given by 

E(?« + Go){qp + Gp)ip-,,2{\q. + G|74e) 



^a/3(q) 



2r(l+p/2)m^Ao 



L G 



G^O 

+ 2r(l + p/2)m, ^ " ■ ^^'^^ <^i+p/2(R?e) 

-2e5„,;3<^p/2(R|e)] . (A. 12) 

Setting p = 1 wc recover Eq. (3.10) of Ref. [28] as we should expect. The case 
p — > is of interest to us. For our vortex lattice we have — 27: pg. Using 

<^o(^)=^e-^ (A.13) 
^^{z)^(l + -^-^e-\ (A.14) 

we see that when letting e — > our dynamical matrix evaluated here agrees 
with the dynamical matrix given in Eq. (26) if we neglect retardation effects. 
We can now go ahead and calculate the 'phonon' or 'magnetophonon' disper- 
sion using 



Here ujc = B/rriy — Tr{ns — nM7)/a^m^ is the 'cyclotron' frequency. Identifying 
the plasma frequency 

as the characteristic frequency wc can evaluate c<Ji(q) and cj2(q) using Eq. (A. 12). 
It turns out that the sums over Rj and G indeed converge rapidly and that 
<^i,2(q) are indeed independent of the value of e. A plot of the spectrum for 
IT'S = nMi (zero 'magnetic' field B) is shown in Fig. 2. Turning on the 'mag- 
netic' field B leads to an avoided crossing of the two modes. This can be seen 
in Fig. 3 where we have chosen lUc = 0.5 Wp. 

Finally we would like to note that in the long- wave length limit the spectrum 
becomes isotopic and we find (with a„ being the distance between nearest 
neighbor vortices) 

2 



^i(q)~4=^M)' (A- 17) 



mp 



uj2{q) ~ ujmp , (A. 18) 
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where Ump — yoUp + is the 'magnetophonon' frequency. With or without 
the 'magnetic' field the shear mode is always linear in q. 



B Simple Debye model and beyond 



It is instructive to evaluate in the Debye approximation where we replace 
cji(q) and a;2(q) by Eqs. (A. 17) and (A. 18). Also, as usual we replace the 
first Brillouin zone by a Debye sphere of the same volume. Using Eq. (32) we 
obtain in this approximation 

^*Ls = (5/2) h/2myUmp , (B.l) 

which for 5 = simplifies to u'^^ = (5/2) h/2myU!p. To this the shear mode 
uJi{q) contributes 80%. Solving for the inertial mass of the vortex we now 
obtain 

1287rV.<s 47rV. " ^ ^ ^ 

More accurately, we can evaluate the integral in Eq. (32) for the exact disper- 
sion relation numerically. Extracting the zero field result we have 

«Ls = 2.5718 • 2^ I{B/myUJp) , (B.3) 

with /(O) = 1. The numerically determined prefactor 2.5718 corresponds to 
the factor 5/2 in the Debye approximation. The small deviation is mainly due 

to the fact that as can be seen in Fig. 2 both cJi(q) and c<j2(q) are on average 
slightly overestimated by Eqs. (A. 17) and (A. 18). Solving for the mass of a 
vortex we find 



(2.5718)^ v/3n^ag , 
- ~^2^^ ^(«r.s^/^) ■ (B.4) 



rms 



If we define h{z) = z^^'^I{z^^'^)/2 then the normahzed function F is related to 
the inverse of I2 by F{x) = a;^/4/^^(a;) and satisfies F{0) = 1. While for the 
case of the exact dispersion relation considered here F{x) has to be calculated 
numerically, Eq. (B.2) suggests a fit of the form F{x) = 1 — CiX^. As can be 
seen in Fig. 4 the quaUty of such a fit turns out to be excellent. 
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